Systematic identification of wheat spike developmental regulators by integrated multi-omics, transcriptional network, GWAS, and genetic analyses.

The spike architecture of wheat plays a crucial role in determining grain number, making it a key trait for optimization in wheat breeding programs. In this study, we used a multi-omic approach to analyze the transcriptome and epigenome profiles of the young spike at eight developmental stages, revealing coordinated changes in chromatin accessibility and H3K27me3 abundance during the flowering transition. We constructed a core transcriptional regulatory network (TRN) that drives wheat spike formation and experimentally validated a multi-layer regulatory module involving TaSPL15, TaAGLG1, and TaFUL2. By integrating the TRN with genome-wide association studies, we identified 227 transcription factors, including 42 with known functions and 185 with unknown functions. Further investigation of 61 novel transcription factors using multiple homozygous mutant lines revealed 36 transcription factors that regulate spike architecture or flowering time, such as TaMYC2-A1, TaMYB30-A1, and TaWRKY37-A1. Of particular interest, TaMYB30-A1, downstream of and repressed by WFZP, was found to regulate fertile spikelet number. Notably, the excellent haplotype of TaMYB30-A1, which contains a C allele at the WFZP binding site, was enriched during wheat breeding improvement in China, leading to improved agronomic traits. Finally, we constructed a free and open access Wheat Spike Multi-Omic Database (http://39.98.48.156:8800/#/). Our study identifies novel and high-confidence regulators and offers an effective strategy for dissecting the genetic basis of wheat spike development, with practical value for wheat breeding.

A large-scale whole-exome sequencing mutant resource for functional genomics in wheat.

Hexaploid wheat (Triticum aestivum), a major staple crop, has a remarkably large genome of ~14.4 Gb (containing 106 913 high-confidence [HC] and 159 840 low-confidence [LC] genes in the Chinese Spring v2.1 reference genome), which poses a major challenge for functional genomics studies. To overcome this hurdle, we performed whole-exome sequencing to generate a nearly saturated wheat mutant database containing 18 025 209 mutations induced by ethyl methanesulfonate (EMS), carbon (C)-ion beams, or γ-ray mutagenesis. This database contains an average of 47.1 mutations per kb in each gene-coding sequence: the potential functional mutations were predicted to cover 96.7% of HC genes and 70.5% of LC genes. Comparative analysis of mutations induced by EMS, γ-rays, or C-ion beam irradiation revealed that γ-ray and C-ion beam mutagenesis induced a more diverse array of variations than EMS, including large-fragment deletions, small insertions/deletions, and various non-synonymous single nucleotide polymorphisms. As a test case, we combined mutation analysis with phenotypic screening and rapidly mapped the candidate gene responsible for the phenotype of a yellow-green leaf mutant to a 2.8-Mb chromosomal region. Furthermore, a proof-of-concept reverse genetics study revealed that mutations in gibberellic acid biosynthesis and signalling genes could be associated with negative impacts on plant height. Finally, we built a publically available database of these mutations with the corresponding germplasm (seed stock) repository to facilitate advanced functional genomics studies in wheat for the broad plant research community.

Boosting wheat functional genomics via an indexed EMS mutant library of KN9204.

A better understanding of wheat functional genomics can improve targeted breeding for better agronomic traits and environmental adaptation. However, the lack of gene-indexed mutants and the low transformation efficiency of wheat limit in-depth gene functional studies and genetic manipulation for breeding. In this study, we created a library for KN9204, a popular wheat variety in northern China, with a reference genome, transcriptome, and epigenome of different tissues, using ethyl methyl sulfonate (EMS) mutagenesis. This library contains a vast developmental diversity of critical tissues and transition stages. Exome capture sequencing of 2090 mutant lines using KN9204 genome-designed probes revealed that 98.79% of coding genes had mutations, and each line had an average of 1383 EMS-type SNPs. We identified new allelic variations for crucial agronomic trait-related genes such as Rht-D1, Q, TaTB1, and WFZP. We tested 100 lines with severe mutations in 80 NAC transcription factors (TFs) under drought and salinity stress and identified 13 lines with altered sensitivity. Further analysis of three lines using transcriptome and chromatin accessibility data revealed hundreds of direct NAC targets with altered transcription patterns under salt or drought stress, including SNAC1, DREB2B, CML16, and ZFP182, factors known to respond to abiotic stress. Thus, we have generated and indexed a KN9204 EMS mutant library that can facilitate functional genomics research and offer resources for genetic manipulation of wheat.

Genome sequencing of Sitopsis species provides insights into their contribution to the B subgenome of bread wheat.

Wheat (Triticum aestivum, BBAADD) is an allohexaploid species that originated from two polyploidization events. The progenitors of the A and D subgenomes have been identified as Triticum urartu and Aegilops tauschii, respectively. Current research suggests that Aegilops speltoides is the closest but not the direct ancestor of the B subgenome. However, whether Ae. speltoides has contributed genomically to the wheat B subgenome and which chromosome regions are conserved between Ae. speltoides and the B subgenome remain unclear. Here, we assembled a high-quality reference genome for Ae. speltoides, resequenced 53 accessions from seven species (Aegilops bicornis, Aegilops longissima, Aegilops searsii, Aegilops sharonensis, Ae. speltoides, Aegilops mutica [syn. Amblyopyrum muticum], and Triticum dicoccoides) and revealed their genomic contributions to the wheat B subgenome. Our results showed that centromeric regions were particularly conserved between Aegilops and Triticum and revealed 0.17 Gb of conserved blocks between Ae. speltoides and the B subgenome. We classified five groups of conserved and non-conserved genes between Aegilops and Triticum, revealing their biological characteristics, differentiation in gene expression patterns, and collinear relationships between Ae. speltoides and the wheat B subgenome. We also identified gene families that expanded in Ae. speltoides during its evolution and 789 genes specific to Ae. speltoides. These genes can serve as genetic resources for improvement of adaptability to biotic and abiotic stress. The newly constructed reference genome and large-scale resequencing data for Sitopsis species will provide a valuable genomic resource for wheat genetic improvement and genomic studies.

Fine mapping of the Hairy glume (Hg) gene in a chromosome variation region at the distal terminus of 1AS.

Trichomes are differentiated epidermal cells and exist on above-ground organs of nearly all land plants with important roles in resistance to a wide range of biotic and abiotic stresses. We attempted to obtain candidate gene (s) for Hairy glume (Hg), responsible for the trichome on wheat glume, by using bulked segregant exome capture sequencing (BSE-Seq), while Hg was only mapped in 0.52-3.26 Mb of 1AS. To further fine map this gene and identify candidate genes in this region, a near isogenic line-derived population consisting of 2,050 F2 lines was generated in the present study. By analyzing this population, Hg was fine mapped into a 0.90 cM region covering a physical distance of ~825.03 Kb encompassing 6 high- and 23 low-confidence genes in the reference genome of Chinese Spring. A presence-absence variation was identified in the fine mapping region through analyses of sequence-tagged sites markers and genome sequences of the hairy glume parent of the near isogenic lines. The results presented here will be useful for further cloning Hg in wheat.

Nutrition Interventions of Herbal Compounds on Cellular Senescence.

When cells undergo large-scale senescence, organ aging ensues, resulting in irreversible organ pathology and organismal aging. The study of senescence in cells provides an important avenue to understand the factors that influence aging and can be used as one of the useful tools for examining age-related human diseases. At present, many herbal compounds have shown effects on delaying cell senescence. This review summarizes the main characteristics and mechanisms of cell senescence, age-related diseases, and the recent progress on the natural products targeting cellular senescence, with the aim of providing insights to aid the clinical management of age-related diseases.

The rice bran peptide KF-8 extends the lifespan and improves the healthspan of Caenorhabditis elegans via skn-1 and daf-16.

With the increased aging of the population, the extension of lifespan and the improvement of healthspan have become important. Our previous studies showed that the rice bran peptide KF-8 exerts an antioxidant effect in cells and mice. In this study, we evaluated the effects of KF-8 on the healthspan and lifespan of Caenorhabditis elegans. We found that KF-8 prolonged the life of nematodes and showed no reproductive toxicity towards nematodes. In addition, KF-8 improved the motility of nematodes and resulted in an extended body length. Using hydrogen peroxide and juglone as stress inducers, we found that KF-8 improved the anti-stress ability of nematodes. In addition, KF-8 upregulated the expressions of skn-1, daf-16 and antioxidant genes. In addition, the life-prolonging effect of KF-8 was lost in skn-1 mutant strains and daf-16 mutant strains, indicating that KF-8 may exert anti-aging effects through skn-1 and daf-16.

Exogenous Bioactive Peptides Have a Potential Therapeutic Role in Delaying Aging in Rodent Models.

In recent years, some exogenous bioactive peptides have been shown to have promising anti-aging effects. These exogenous peptides may have a mechanism similar to endogenous peptides, and some can even regulate the release of endogenous active peptides and play a synergistic role with endogenous active peptides. Most aging studies use rodents that are easy to maintain in the laboratory and have relatively homogenous genotypes. Moreover, many of the anti-aging studies using bioactive peptides in rodent models only focus on the activity of single endogenous or exogenous active peptides, while the regulatory effects of exogenous active peptides on endogenous active peptides remain largely under-investigated. Furthermore, the anti-aging activity studies only focus on the effects of these bioactive peptides in individual organs or systems. However, the pathological changes of one organ can usually lead to multi-organ complications. Some anti-aging bioactive peptides could be used for rescuing the multi-organ damage associated with aging. In this paper, we review recent reports on the anti-aging effects of bioactive peptides in rodents and summarize the mechanism of action for these peptides, as well as discuss the regulation of exogenous active peptides on endogenous active peptides.

High animal protein diet and gut microbiota in human health.

The role of the intestinal flora in health and disease has become a research hotspot. Compared with carbohydrates and fats, proteins are metabolized primarily by microbial fermentation in the intestine. The production of protein fermentation products and metabolites depends on the composition, diversity, and metabolism of the gut microbiota. Several protein fermentation products, including indoles, phenols, polyamines, hydrogen sulfide (H2S), amines, and carnitine, are toxic. This study analyzes the relationship between high-protein diets (HPDs), the intestinal microbiota, and human health and disease. Long-term HPDs increase the risk of intestinal diseases, type 2 diabetes (T2DM), obesity, central nervous system (CNS) diseases, and cardiovascular diseases (CVD) by producing toxic metabolites in the colon, including amines, H2S, and ammonia. Short-term HPDs have little effect on the metabolism of healthy individuals under 65 years old. However, meeting the protein requirements of individuals over 65 years old using HPDs is more challenging. The adverse effects of HPDs on athletes are minimal. Natural compounds (plant extracts, whose main constituents are polysaccharides and polyphenols), prebiotics, probiotics, and regular physical exercise improve gut dysbiosis and reduce disease risk.

Design and Analysis of a Radio-Frequency Moisture Sensor for Grain Based on the Difference Method.

Grain moisture is one of the key indexes of grain quality, and acquiring an accurate moisture value is critical for grain storage security. However, the sensors used in the traditional methods for testing grain moisture are based on capacitance, microwave, or radio-frequency methods and still exhibit low accuracy and instability because they are susceptible to the temperature, moisture, and micro gas flow of the air in the granary. In this study, we employed a new design for a radio-frequency moisture sensor for grain. The structure of the sensor is based on the difference method and consists of two parallel probe units. These units are at different distances to the tested grain, resulting in different sensitivities in the moisture measurements. Through a phase difference operation on the test signals, the disturbance variable was reduced. The specific size of the two parallel probes was confirmed by calculation and simulation using High Frequency Structure Simulator (HFSS) software. The simulated and measured parameters of a prototype sensor agreed well. The linear relationship yielded a correlation coefficient of 0.9904, and the average error of the moisture testing was within ±0.3% under the conditions where the VSWR (voltage standing wave ratio) value and return losses were 1.5896 and -20 dB, respectively, at a measured central frequency of 100 MHz. The results indicate that the performance of the sensor was excellent.

A pig BodyMap transcriptome reveals diverse tissue physiologies and evolutionary dynamics of transcription.

A comprehensive transcriptomic survey of pigs can provide a mechanistic understanding of tissue specialization processes underlying economically valuable traits and accelerate their use as a biomedical model. Here we characterize four transcript types (lncRNAs, TUCPs, miRNAs, and circRNAs) and protein-coding genes in 31 adult pig tissues and two cell lines. We uncover the transcriptomic variability among 47 skeletal muscles, and six adipose depots linked to their different origins, metabolism, cell composition, physical activity, and mitochondrial pathways. We perform comparative analysis of the transcriptomes of seven tissues from pigs and nine other vertebrates to reveal that evolutionary divergence in transcription potentially contributes to lineage-specific biology. Long-range promoter-enhancer interaction analysis in subcutaneous adipose tissues across species suggests evolutionarily stable transcription patterns likely attributable to redundant enhancers buffering gene expression patterns against perturbations, thereby conferring robustness during speciation. This study can facilitate adoption of the pig as a biomedical model for human biology and disease and uncovers the molecular bases of valuable traits.

Myb10-D confers PHS-3D resistance to pre-harvest sprouting by regulating NCED in ABA biosynthesis pathway of wheat.

Pre-harvest sprouting (PHS), the germination of grain before harvest, is a serious problem resulting in wheat yield and quality losses. Here, we mapped the PHS resistance gene PHS-3D from synthetic hexaploid wheat to a 2.4 Mb presence-absence variation (PAV) region and found that its resistance effect was attributed to the pleiotropic Myb10-D by integrated omics and functional analyses. Three haplotypes were detected in this PAV region among 262 worldwide wheat lines and 16 Aegilops tauschii, and the germination percentages of wheat lines containing Myb10-D was approximately 40% lower than that of the other lines. Transcriptome and metabolome profiling indicated that Myb10-D affected the transcription of genes in both the flavonoid and abscisic acid (ABA) biosynthesis pathways, which resulted in increases in flavonoids and ABA in transgenic wheat lines. Myb10-D activates 9-cis-epoxycarotenoid dioxygenase (NCED) by biding the secondary wall MYB-responsive element (SMRE) to promote ABA biosynthesis in early wheat seed development stages. We revealed that the newly discovered function of Myb10-D confers PHS resistance by enhancing ABA biosynthesis to delay germination in wheat. The PAV harboring Myb10-D associated with grain color and PHS will be useful for understanding and selecting white grained PHS resistant wheat cultivars.

Combining a New Exome Capture Panel With an Effective varBScore Algorithm Accelerates BSA-Based Gene Cloning in Wheat.

The discovery of functional genes underlying agronomic traits is of great importance for wheat improvement. Here we designed a new wheat exome capture probe panel based on IWGSC RefSeq v1.0 genome sequence information and developed an effective algorithm, varBScore, that can sufficiently reduce the background noise in gene mapping and identification. An effective method, termed bulked segregant exome capture sequencing (BSE-Seq) for identifying causal mutations or candidate genes was established by combining the use of a newly designed wheat exome capture panel, sequencing of bulked segregant pools from segregating populations, and the robust algorithm varBScore. We evaluated the effectiveness of varBScore on SNP calling using the published dataset for mapping and cloning the yellow rust resistance gene Yr7 in wheat. Furthermore, using BSE-Seq, we rapidly identified a wheat yellow leaf mutant gene, ygl1, in an ethyl methanesulfonate (EMS) mutant population and found that a single mutation of G to A at 921 position in the wild type YGL1 gene encoding magnesium-chelatase subunit chlI caused the leaf yellowing phenotype. We further showed that mutation of YGL1 through CRISPR/Cas9 gene editing led to a yellow phenotype on the leaves of transgenic wheat, indicating that ygl1 is the correct causal gene responsible for the mutant phenotype. In summary, our approach is highly efficient for discovering causal mutations and gene cloning in wheat.

Genome-wide DNA methylation profile of prepubertal porcine testis.

The biological structure and function of the mammalian testis undergo important developmental changes during prepuberty and DNA methylation is dynamically regulated during testis development. In this study, we generated the first genome-wide DNA methylation profile of prepubertal porcine testis using methyl-DNA immunoprecipitation (MeDIP) combined with high-throughput sequencing (MeDIP-seq). Over 190 million high-quality reads were generated, containing 43642 CpG islands. There was an overall downtrend of methylation during development, which was clear in promoter regions but less so in gene-body regions. We also identified thousands of differentially methylated regions (DMRs) among the three prepubertal time points (1 month, T1; 2 months, T2; 3 months, T3), the majority of which showed decreasing methylation levels over time. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that many genes in the DMRs were linked with cell proliferation and some important pathways in porcine testis development. Our data suggest that DNA methylation plays an important role in prepubertal development of porcine testis, with an obvious downtrend of methylation levels from T1 to T3. Overall, our study provides a foundation for future studies and gives new insights into mammalian testis development.

Uncovering the dispersion history, adaptive evolution and selection of wheat in China.

Wheat was introduced to China approximately 4500 years ago, where it adapted over a span of time to various environments in agro-ecological growing zones. We investigated 717 Chinese and 14 Iranian/Turkish geographically diverse, locally adapted wheat landraces with 27 933 DArTseq (for 717 landraces) and 312 831 Wheat660K (for a subset of 285 landraces) markers. This study highlights the adaptive evolutionary history of wheat cultivation in China. Environmental stresses and independent selection efforts have resulted in considerable genome-wide divergence at the population level in Chinese wheat landraces. In total, 148 regions of the wheat genome show signs of selection in at least one geographic area. Our data show adaptive events across geographic areas, from the xeric northwest to the mesic south, along and among homoeologous chromosomes, with fewer variations in the D genome than in the A and B genomes. Multiple variations in interdependent functional genes such as regulatory and metabolic genes controlling germination and flowering time were characterized, showing clear allelic frequency changes corresponding to the dispersion of wheat in China. Population structure and selection data reveal that Chinese wheat spread from the northwestern Caspian Sea region to South China, adapting during its agricultural trajectory to increasingly mesic and warm climatic areas.

Identification and positional distribution analysis of transcription factor binding sites for genes from the wheat fl-cDNA sequences.

The binding sites of transcription factors (TFs) in upstream DNA regions are called transcription factor binding sites (TFBSs). TFBSs are important elements for regulating gene expression. To date, there have been few studies on the profiles of TFBSs in plants. In total, 4,873 sequences with 5' upstream regions from 8530 wheat fl-cDNA sequences were used to predict TFBSs. We found 4572 TFBSs for the MADS TF family, which was twice as many as for bHLH (1951), B3 (1951), HB superfamily (1914), ERF (1820), and AP2/ERF (1725) TFs, and was approximately four times higher than the remaining TFBS types. The percentage of TFBSs and TF members showed a distinct distribution in different tissues. Overall, the distribution of TFBSs in the upstream regions of wheat fl-cDNA sequences had significant difference. Meanwhile, high frequencies of some types of TFBSs were found in specific regions in the upstream sequences. Both TFs and fl-cDNA with TFBSs predicted in the same tissues exhibited specific distribution preferences for regulating gene expression. The tissue-specific analysis of TFs and fl-cDNA with TFBSs provides useful information for functional research, and can be used to identify relationships between tissue-specific TFs and fl-cDNA with TFBSs. Moreover, the positional distribution of TFBSs indicates that some types of wheat TFBS have different positional distribution preferences in the upstream regions of genes.

Genome-Wide Association Study for Pre-harvest Sprouting Resistance in a Large Germplasm Collection of Chinese Wheat Landraces.

Pre-harvest sprouting (PHS) is mainly caused by the breaking of seed dormancy in high rainfall regions, which leads to huge economic losses in wheat. In this study, we evaluated 717 Chinese wheat landraces for PHS resistance and carried out genome-wide association studies (GWAS) using to 9,740 DArT-seq and 178,803 SNP markers. Landraces were grown across six environments in China and germination testing of harvest-ripe grain was used to calculate the germination rate (GR) for each accession at each site. GR was highly correlated across all environments. A large number of landraces (194) displayed high levels of PHS resistance (i.e., mean GR < 0.20), which included nine white-grained accessions. Overall, white-grained accessions displayed a significantly higher mean GR (42.7-79.6%) compared to red-grained accessions (19.1-56.0%) across the six environments. Landraces from mesic growing zones in southern China showed higher levels of PHS resistance than those sourced from xeric areas in northern and north-western China. Three main quantitative trait loci (QTL) were detected by GWAS: one on 5D that appeared to be novel and two co-located with the grain color transcription factor Tamyb10 on 3A and 3D. An additional 32 grain color related QTL (GCR-QTL) were detected when the set of red-grained landraces were analyzed separately. GCR-QTL occurred at high frequencies in the red-grained accessions and a strong correlation was observed between the number of GCR-QTL and GR (R2 = 0.62). These additional factors could be critical for maintaining high levels of PHS resistance and represent targets for introgression into white-grained wheat cultivars. Further, investigation of the origin of haplotypes associated with the three main QTL revealed that favorable haplotypes for PHS resistance were more common in accessions from higher rainfall zones in China. Thus, a combination of natural and artificial selection likely resulted in landraces incorporating PHS resistance in China.

SNP Mining in Functional Genes from Nonmodel Species by Next-Generation Sequencing: A Case of Flowering, Pre-Harvest Sprouting, and Dehydration Resistant Genes in Wheat.

As plenty of nonmodel plants are without genomic sequences, the combination of molecular technologies and the next generation sequencing (NGS) platform has led to a new approach to study the genetic variations of these plants. Software GATK, SOAPsnp, samtools, and others are often used to deal with the NGS data. In this study, BLAST was applied to call SNPs from 16 mixed functional gene's sequence data of polyploidy wheat. In total 1.2 million reads were obtained with the average of 7500 reads per genes. To get accurate information, 390,992 pair reads were successfully assembled before aligning to those functional genes. Standalone BLAST tools were used to map assembled sequence to functional genes, respectively. Polynomial fitting was applied to find the suitable minor allele frequency (MAF) threshold at 6% for assembled reads of each functional gene. SNPs accuracy form assembled reads, pretrimmed reads, and original reads were compared, which declared that SNPs mined from the assembled reads were more reliable than others. It was also demonstrated that mixed samples' NGS sequences and then analysis by BLAST were an effective, low-cost, and accurate way to mine SNPs for nonmodel species. Assembled reads and polynomial fitting threshold were recommended for more accurate SNPs target.

Genome-wide characterization of developmental stage- and tissue-specific transcription factors in wheat.

BACKGROUND: Wheat (Triticum aestivum) is one of the most important cereal crops, providing food for humans and feed for other animals. However, its productivity is challenged by various biotic and abiotic stresses such as fungal diseases, insects, drought, salinity, and cold. Transcription factors (TFs) regulate gene expression in different tissues and at various developmental stages in plants and animals, and they can be identified and classified into families according to their structural and specialized DNA-binding domains (DBDs). Transcription factors are important regulatory components of the genome, and are the main targets for engineering stress tolerance. RESULTS: In total, 2407 putative TFs were identified from wheat expressed sequence tags, and then classified into 63 families by using Hmm searches against hidden Markov model (HMM) profiles. In this study, 2407 TFs represented approximately 2.22% of all genes in the wheat genome, a smaller proportion than those reported for other cereals in PlantTFDB V3.0 (3.33%-5.86%) and PlnTFDB (4.30%-6.46%). We assembled information from the various databases for individual TFs, including annotations and details of their developmental stage- and tissue-specific expression patterns. Based on this information, we identified 1257 developmental stage-specific TFs and 1104 tissue-specific TFs, accounting for 52.22% and 45.87% of the 2407 wheat TFs, respectively. We identified 338, 269, 262, 175, 49, and 18 tissue-specific TFs in the flower, seed, root, leaf, stem, and crown, respectively. There were 100, 6, 342, 141, 390, and 278 TFs specifically expressed at the dormant seed, germinating seed, reproductive, ripening, seedling, and vegetative stages, respectively. We constructed a comprehensive database of wheat TFs, designated as WheatTFDB ( http://xms.sicau.edu.cn/wheatTFDB/ ). CONCLUSIONS: Approximately 2.22% (2407 genes) of all genes in the wheat genome were identified as TFs, and were clustered into 63 TF families. We identified 1257 developmental stage-specific TFs and 1104 tissue-specific TFs, based on information about their developmental- and tissue-specific expression patterns obtained from publicly available gene expression databases. The 2407 wheat TFs and their annotations are summarized in our database, WheatTFDB. These data will be useful identifying target TFs involved in the stress response at a particular stage of development.